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planar limit could display a phase transition in that limit as nonlinear effects become 
dominating over linear ones. 
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1. Introduction 



1.1 Generalities 

Recent numerical work provides evidence tliat Wilson loops in SU{N) gauge theory in two, 
three and four dimensions exhibit an infinite- phase transition as they are dilated from 
a small size to a large one. In the course of this dilation the eigenvalue distribution of 
the untraced Wilson loop unitary matrix expands from a small arc on the unit circle to 
encompassing the entire unit circle [1, 2]. The universality class of this transition is that 
of a random multiplicative ensemble of unitary matrices. The transition was discovered by 
Durhuus and Olesen [3] (DO) when they solved the Migdal-Makeenko [4] loop equations in 
two dimensional planar QCD. The associated multiplicative random matrix ensemble [5] 
can be axiomatized in the language of noncommutative probability [6] . It provides a gen- 
eralization of the familiar law of large numbers. The essential feature making a difference 
is that one case is commutative and the other is not. Various recent insights into the 
DO transition [7, 8, 9] indicate that an even deeper understanding of the transition might 
emerge. The large- transition does not imply confinement, but if confinement occurs 
the large- transition is unavoidable if parallel transport round small loops is close to the 
identity. 

In this paper we relax the unitarity constraint. We shall focus on a simple multiplicative 
random complex matrix model introduced in [10], where it was shown that the model has 
an infinite- phase transition. The motivating physics application for this study is to a 
more general gauge theory, obeying extra symmetries, which make complex matrix valued 
Wilson loop operators natural observables. If the situation outlined above for ordinary 
gauge theories, where the matrix of the Wilson loop operator is unitary, generalizes to the 
complex case, the multiplicative random complex matrix model might capture the universal 
features of large- A^ transitions occurring in these more elaborate models. 

We also wish to point out that complex matrix transitions may also be relevant to 
ordinary gauge theories, in dimensions higher than two, for technical reasons. Ultraviolet 
divergences of the bare Wilson loop matrix may be eliminated by a regularization prescrip- 
tion that makes the Wilson loop operator non-unitary. For concreteness, let us assume we 
are in four dimensions, dealing with planar QCD. One example of a regularization is to 
associate with a curve C that has a marked point x on it the operator 

liA-dy+<f{y)\dy\] 

W{C,x)=Pec . (1.1) 

The ordered line integral starts at x and follows the oriented curve C until it gets back to 
X. We take the extra scalar field $ = $"1" transforming as an adjoint under the gauge group, 
with a mass much heavier than the QCD scale Aqcd- A = is the usual gauge field, and 
y^(s) describes the curve C. By adjusting the normalization of its contribution could 
be made to cancel out the linear perimeter divergence associated with W but otherwise 
have little impact on smooth loops larger than the QCD scale on account of its large mass. 
For this to work $ must enter the exponent without a \/— T prefactor. If the mass M of 
$ is very large, one could extract the pure QCD string tension from loops of area A with 
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1 ^ AqcdV^ ^ M/Aqcd- The rcgularization would make W a finite operator, but one 
cannot associated it with a unitary matrix, and its spectrum would be spread somewhat 
in the complex plane, defining a finite surface eigenvalue density rather than a finite linear 
eigenvalue density. 

1.2 Outline of the paper 

We start by presenting a set of symmetries we require the complex Wilson loop matrices 
to obey. These symmetries are inherited from the Euclidean Gauge Theory producing 
the Wilson loops. We argue that these symmetries, in conjunction with the assumption 
that the Wilson loop W has a perturbative weak-coupling regime and a non-perturbative 
"disordered" regime, point to a largc-iV phase transition in the spectrum of W. The 
support of the eigenvalues of W undergoes a topological change at the transition point, 
and this indicates that there is something universal about the transition. We shall refer 
to this hypothetical universality as large-iV universality, to distinguish it from ordinary 
critical phenomena universality, which governs the structure of the pertinent Euclidean 
Gauge Theory at any finite N and, by assumption, extends smoothly to infinite N. 

Armed with the universality assumption we then make a guess for the simplest possible 
random matrix model which would be in the same large- universality class as the above 
Wilson loops. We proceed by discussing the general properties of the model. Some time 
is devoted to a technical point: Simplifications occur when one drops the det W = 1 
constraint, but as in usual planar QCD, dropping the constraint has no impact on the 
infinite-AT phase structure. In the context of the model we can be more specific, without 
actually solving the model, about the shape of the support of the spectrum of W at infinite 
A^. We find that it is constrained to an annulus in the complex plane whose internal and 
external radius are reciprocal. As the coupling or loop size of W change the spectrum 
evolves from a simply connected small blob centered at 2; = 1 to a multiply connected 
region. The annulus the spectrum is confined to also expands. This picture mirrors the 
picture we presented in the Euclidean Gauge Theory case described earlier. 

We next proceed to a more detailed analysis of the model using the average of the mod- 
ulus square of the characteristic polynomial of W. This observable is shown quantitatively 
to produce a spectrum with the properties that were anticipated. We use analytical and 
also some numerical tools. The main analytical tool is a representation of the observable 
in terms of an integral over Grassmann variables with a local action in an internal space; 
the Grassmann variables are akin to mathematical quarks. After the introduction of the 
Grassmann variables the matrix averaging of the model can be done explicitly and, even- 
tually, the entire dependence on A^ becomes explicit. This sets the stage for a saddle-point 
analysis at infinite A^. We perform the analysis only to the extent that it gives the phase 
structure. Global stability questions are dealt with by numerical tests and not by purely 
analytical methods. 

To get some feeling for the universality of the model we proceed with a slight general- 
ization. This generalization has an extra parameter which allows an interpolation between 
the complex multiplicative matrix model and the unitary multiplicative matrix model that 
has been extensively studied in past work. In that, the generalized model provides further 
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support to the view that in some sense the large- transition here has a direct relation- 
ship to the large-A'' transition found in ordinary gauge theories with unitary Wilson loop 
matrices. It is seen that, similarly to the original complex matrix model and to the unitary 
matrix model, the inviscid Burgers equation plays a central role also in the generalized 
model. 

The last part of the paper is only partially successful. Although the infinite- A'^ phase 
structure justifies the guess that there is a large-A^ universality class associated with it, 
to make this concrete one needs to go to sub-leading terms in the large- iV expansion and 
identify the (hopefully few) relevant variables. We note that for finite N the matrix models 
get mapped into matrix models consisting of products of 2 x 2 matrices. More precisely, 
the average of the modulus square of the characteristic polynomial of the N x N complex 
matrix W can be exactly represented by the solution of an associated multiplicative matrix 
model where the matrices are only 2x2 and the dependence on N is explicit. However, 
this still leaves too many variables (albeit a finite, A^-independent, number), impeding 
an explicit analysis of the approach to the large- A/" limit. Suspecting that the number of 
variables can be further reduced by dropping sub-leading corrections in large N we simplify 
the 2x2 model further focusing on some special cases. We finally present a case where 
we end up with only two real variables and show how that model could be exactly solved. 
However, the exact solution is in the form of an infinite series, and the study of the large- A'^ 
limit still presents difficulties. 

We end up being forced to leave further work on the large- A?^ universality class to 
the future, but feel that we have made substantial progress here to eventually achieve a 
complete understanding of the large- AT universality class. 

The paper ends with a brief summary. 

2. General properties of our random matrix products 

The complex Wilson loop matrix is denoted by W. We assume that the Euclidean Field 
Theory which defines W provides a probability distribution for W, P{W)(P^^~'^W , with 
some natural properties: 

• dctW = 1, 

• P{W) = P{W-^) = P{W*), 

• P{W) = PiUWU'f) for U G U{N). 

A construction of W in terms of traceless double indexed fields and discrete symmetries 
like parity and charge conjugation can assure the first two properties. Gauge invariance 
implies the third property. 

To study the spectral properties of W^, we define 



Qiz,z*) = {\detiz-W)\''), 



(2.1) 



where (...) denotes averaging with respect to P. The general properties imply 



Q{z,z*) = Q{z\z) 



2N 



Q{l/z*,l/z) 



= \z 



2N 



Q{l/z,l/z*). 



(2.2) 



- 5 - 



P and consequentially Q are assumed to depend on a coupling A > (A can be a running 
coupling, depending on the size of the loop). When — > oo A will stay finite, scaling with 
N to produce the 't Hooft topological classification of diagrams. We assume that we have 
employed a regularization that respects the above symmetries, admits a standard large- 
expansion and is in general benign in the sense that a large- transition will survive the 
continuum limit. 

Generically, N x N random matrices have regions in the spectral plane where the 
eigenvalue density is exponentially suppressed as A'' ^ oo; hence, at A'" = oo the eigenvalue 
density vanishes in these regions. The eigenvalue density at infinite A^ cannot vanish 
everywhere, so the large- A^ limit induces some lack of smoothness in the eigenvalue density. 
Typically, the eigenvalue density is guaranteed to be non-zero somewhere in the plane 
(rather than disappearing at infinity) because for a small loop, W is close to the identity 
matrix. In many theories one can replace the phrase "small loop" by "small A" . However, 
the condition det W = 1 ensures that is not an eigenvalue. 

We therefore assume that in the N = oo limit, for any A 7^ the probability of W 
having eigenvalues within some small finite circle around 2; = is zero. The radius of 
the circle increases to unity when A — > when all the eigenvalues of W are forced into a 
shrinking region around z = 1. Typically, this is reflected in Q{z, z*) having a holomorphic 
factorized form for l^l < p(A) < 1, 



By the inversion symmetry eigenvalues are also excluded from around complex infinity, 
so the complex plane can be thought of as a two-dimensional sphere with the north and 
south poles excised. At infinite A^ the eigenvalues make up a connected region containing 
z = 1 for any A. For A <C 1 this region is very small and does not wrap around the doubly 
punctured sphere. 

When A 00, the dynamics of the particular model become important. We are 
interested here in the case where strong coupling induces strong disorder in the Wilson 
loop spectrum. This is the case in confining theories, but the evolution to sufficiently 
strong disorder causing a large- A^ transition is not a compelling reason for the onset of 
confinement in more general theories. Strong disorder would imply that as the coupling 
becomes stronger, P{W) becomes less restrictive, and the set of eigenvalues of W spreads 
widely. In this case it makes sense to assume the spectrum to completely surround the 
origin z = 0. As a result the simply connected domain where the eigenvalues reside at 
small couplings becomes topologically nontrivial on the doubly punctured sphere, becoming 
multiply connected as a result of the punctures. In principle, more complicated topology 
changes could happen, but, under some conditions, extra restrictions similar to the result 
of [11] might apply, leaving us with only the simplest option just described. Intuitively, 
this is the generic way in which eigenvalues would spread out as disorder increases in a 
model obeying the general symmetries described earlier. 

As already mentioned, z = 1 will typically be in the domain of eigenvalues. Thus, 
the unit circle intersects the set of possible eigenvalues and we could look for a signal of 
the transition on |2;| = 1. This signal would be the entrance of the point z = —1 into the 




(2.3) 
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domain of eigenvalues as A is increased from through the transition point. The points 
z = ±1 are special because they are fixed points of the inversion symmetry, z = 1 is in the 
domain of eigenvalues for any A, and the transition occurs when z = —1 also joins. This 
description makes the similarity with the unitary matrix case clear. 

3. General considerations about a basic random complex matrix model 

We now set up a simple random matrix model for the complex Wilson loop matrix. Basi- 
cally we replace the true P(W) by a much simpler one. The model is almost identical to 
that of [10]; the minor difference is irrelevant, but shall nevertheless be touched upon later. 
The integration measure over complex numbers z = x + iy is defined as 

dfi{z) = 2^'^^ ~ dxdy . (3.1) 

3.1 SL{N,C) case 

Consider the space of traceless N x N complex matrices C and define a normalized prob- 
ability density over it, 

P{C)d^i{C) = e-^''^^^TrS{tiC) ff -df,{Cij) , (3.2) 

J. J. yj- 

where the complex delta function is defined as 6{z) = 6{x)5{y). For any complex matrices 
A and B we have 

J P{C)dfi{C)e''^'^+'''''^ = e^trBtA-^trAtrBt ^33^ 

Define now a sequence of n i.i.d. matrices Mj, j = 1, . . . ,n, 

Mj = e'^^ , (3.4) 

where Cj is distributed by P{Cj) and e > is a small number. The delta function in the 
probability density (3.2) ensures det(Mj) = 1 for all j = 1, . . . ,n. The distributions of the 
Cj's are invariant individually under Cj — > C*, —Cj and Cj — > UjCjUj, Uj G U{N). 
Define 

n 

Wn = Ml Ms ■■■Mn = '[lMj. (3.5) 

i=i 

Mj,M*,M~^ and WMjU, U £ U{N) are equally probable. Also, any two permuted 
sequences of M^'s are equally probable. The distribution of Wn has the properties listed 
in the previous section. 

We are interested in the limit n — > oo, e — > with t = e^n held fixed at a non-negative 
value. ^ In that limit the product matrix W will be a finite matrix, and we are interested 
in properties of its distribution as a function of t. 

^£ could have been complex, but C — > e**C can be used to make e > 0. 
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3.2 GL{N,C) case 

The matrices W were strictly restricted to have unit determinant. Imposing the hnear 
restriction trC = forces Wn G SL{N,C) for any sequence Mj. However, Httle is lost 
by relaxing the determinant restriction. As we shall see later, this has no effect on the 
boundary of the region of non- vanishing eigenvalue density in the infinite- limit. 

Without the restriction on the determinant, we would define a probability density 
(instead of (3.2)) by 



P{C)d^^{C)=e-''''^'^ TT -dfi{Qj), (3.6) 

J. J. Yj- 

l<ij<iV 



which would then lead to 

J P(C)d/x(C)e*'-^'^+*'-^'^ = 6^*=-^'^ . (3.7) 

In this case, we are left with only the first term in the exponent of (3.3). 

3.3 Fokker-Planck equation and why det W = 1 does not matter at infinite N 

In this section we drop the restriction tr C = and explain why this does not affect the 
N = oo limit. The joint probability distribution of the entries of the matrix Wn in the 
limit t fixed, for finite N, is PNiWa/s] t)d^{W) , where the Wap are iV^ complex numbers 
and dfi{W) is some conveniently chosen measure which is independent of t. 

We have a Markov chain [12, 13], and P/v will satisfy a partial differential equation of 
the form [14] 

^ = ^nPn , (3.8) 

where J7jv is linear partial differential operator of at most second order acting on the 2A^ 
real variables defining the A^ complex numbers Waf3- depends explicitly on A" but has 
no explicit t-dependence. 

r^TV is determined by the terms of order e and in M„ in the recursion 

Wn = Wn-lMn . (3.9) 

Higher-order terms do not matter. (3.8) is an equation of the Fokker-Planck (FP) type. 
P/v is determined by the initial condition at t = 0, which in our case is a 5- function with 
respect to d^iW) concentrated at W given by the unit matrix. 

The FP equation for P/v is derived by first expressing the step-n probability density 
P^\w)d^x{W) in terms of P^~^\w')d^{W'), where is kept fixed. The computation 
is based on the linear recursion (3.9), 

W = W'M , (3.10) 

where M = e^*^ with C distributed according to P{C)dn{C). One expresses P^\w)dii{W) 
in terms of W' going to order e^. M acts linearly on the rows of W' ^ and therefore 
the Jacobian is given by the product of the Jacobians per row, which is |detMp^ = 
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g2eRetrC_ 'pj^jg takes care of the change in the measure. The expansion in W' around W 
produces first-order derivatives of at order e and and second-order derivatives 

at order e^. After averaging over the matrix C, the measure terms of order give terms 
proportional to pj^ ^\ while measure terms of order e can combine with terms of order e 
from the expansion of P^ giving first-order derivative terms in P^ ^\ All second-order 
derivative terms in P^ come from the expansion of and not from the measure. 

This discussion simplifies if one chooses a measure term that is invariant under the 
recursion. This is possible in many cases when the evolution is on a group manifold. In our 
case it is convenient to parametrize W S GL{N,C) as W = wW with = det(VF) and 
W G SL{N,C). Correspondingly we factorize the measure, diJ,{W) = d^{w)dfi{W), where 
diJ,{W) is right-invariant on SL{N,C). P^\w) also factorizes in these variables. This 
follows by induction in n, since if P^ is factorized so is P^^ , and the initial condition 
also factorizes. Hence 

PPiW) = Pj^^\w)Pi,^^\w) . (3.11) 

P^^^ is very easy to compute, as it comes from an abelian ensemble. Because of the choice 
for the measure term and because of the invariance of M under conjugation by SU{N) 
elements, the FP equation is invariant under W WU with U £ SU{N). 

One can now restore the relation between n and e and take e ^ 0, n ^ oo at fixed t, 
leading to 

pP{w, W) ^ pn{w; t)PN{W; t) . (3.12) 

If W were a product of SU{N) matrices, because of the invariance under right multipli- 
cation by elements of SU{N), the form of the operator ^Ijy would be uniquely fixed up 
to an overall constant to being the Laplacian on the SU{N) group manifold, and the FP 
equation would become the heat-kernel equation on SU{N), leaving only an overall scale 
to be determined. 

Parametrizing det W = e""^™ with n, f € M, it turns out that both u and v are normally 
distributed and their distributions are A^-in dependent. Therefore, replacing by in 
averages of characteristic polynomials (cf. Sec. 4) only produces an unimportant prefactor 
in the large-A limit. As a result, we see that we can use GL{N,C) instead of SL{N,C) 
without affecting the large- A^ limit. At subleading order in 1/A there are differences, but 
they are easily determined from pi\i{'w;t). So, little is lost by working with GL{N,C) 
instead of SL{N, C). 

3.4 The limit n — oc, e — > Q, t = n^e fixed and the subsequent N — ^ oo Hmit 

Our objective is to determine the region in the z- plane populated by eigenvalues of Wn- As 
a first step, we would like to find some bounds delimiting this region. The region will have 
a sharp boundary after the A — > oo limit is taken. Even without the restriction det W = 1, 
at large A^, exact inversion symmetry gets restored. We expect that after all the limits are 
taken we shall have 

e-^W < |z| < e^W, z G spectrum(W„) . (3.13) 
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The function 'y(t) is positive for all t > 0. We expect 7(t) to be finite for all finite t and to 
increase monotonically with t, because of the associated increase in disorder. By definition, 
at any fixed t > 0, 7(t) is the smallest positive number for which (3.13) holds almost surely, 
i.e., with probability 1. The annulus keeps the spectrum of Wn away from the origin and 
infinity for any finite t. The major structural change that can happen as t is increased from 
is that the spectrum wraps round the annulus. For small t the spectrum is a small blob 
round z = 1; the inversion and reflection invariance give it a kidney shaped appearance. 
As t increases the blob has a larger annulus available and expands into it, until eventually 
it reaches around it, at some finite critical t. 

To get some feeling for why there is a j{t) at all and how it behaves we start from some 
small e, large n and large A^, without committing at the moment to any special relations 
between these numbers, although we really are interested in the situation n ~ 
and, although ^ 1, we want N. 

3.4.1 t^oc 

We first fix some A ^ 1, take some fixed e'^ <^ 1 and study what happens as n — > oo. 
In terms of our true interest this means we are are trying to understand the asymptotic 
behavior of 7(i) for t = ne^ going to infinity. 

If we take n — > oo at fixed > the classical results of Fiirstenberg [15, 16] and a 
theorem by Oseledec [17] apply. Some standard texts on the topic of random multiplicative 
matrix ensembles are [18, 19]. Our discussion is generally based on [15, 16, 18, 19] and 
specifically on [20, 21], but we do not aim here for mathematical rigor. 

We start by defining the norm of a vector v G by 

N 

(3.14) 

a=l 

and the matrix norm by^ 

llWnll = sup \\vWn\\ , (3.15) 
lbll=l 

identifying it with the square root of the largest eigenvalue of WnWn. The results about 
random matrix products mentioned above apply to Wn- For a fixed Wn we have 

llWnll > e^W . (3.16) 

The inequality may be sharp in the above equation because z, the eigenvalue of Wn that 
has maximum absolute value and defines ■j{t), could be associated to an eigenvector that is 
very different from the maximum eigenvector of WnWn for all asymptotic times t. However, 
we expect that our case is generic enough for a conjecture of [22] to apply, which would 
allow us to replace the inequality sign above by an asymptotic equality at infinite t. For a 
discussion, see [19], page 21. If this is assumed, (3.16) can be replaced by an equality. 

^In Dirac notation, vW„ — > {i;|M^„. 
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This assumption is non-trivial: For example, in the Ginibre ensemble [23], the left- 
hand side of (3.16) equals twice the right-hand side. However, in this case Wn is not given 
by a product, but is just a complex matrix C distributed according to exp[— A'"tr(C^^C)], 
and no noncommutative matrix products are involved.'^ This is very different from Wn 
for large t, but is intuitively close to the situation for small t, the case addressed in the 
next subsection. There, our estimate for -y{t) will be direct, without involving the norm 
||Wn||- In any case, that (3.16) becomes an equality as t — > oo will be confirmed by both 
the analytical and numerical results presented later in the paper. 

To study II Wnll we need to know what happens to vWn for an arbitrary v with ||?; || / 0. 
We define Vi, i = 1, . . . ,n hy 

i 

Vi = vllMj (3.17) 

i=i 

and vq = V. We are only interested in the ray specified by v. Let 

||f II 

Because of the invariance under conjugation by U{N) elements, the distribution of 
induced by that of Mi is independent of v. We now write 

n 11 11 

logll^nll -log||v|| = logf^ (3.19) 

1 \\Vm~l\\ 
m=l 

This is a trick used in [20] . 

The terms in the sum on the RHS of (3.19) are i.i.d. real numbers for any fixed values 
of n, e^, by the same argument as below Eq. (3.18). Therefore, we can calculate the 
probability distribution of the LHS by calculating the characteristic function F{k) (see 
Eq. (3.22) below) of one of the terms on the RHS, taking the power n, and taking the 
inverse Fourier transform of that. 

We reproduce the equation describing the source of randomness, ignoring the zero 
trace condition, 

M = e'^, P{C)=Me-^^'^^^ . (3.20) 
The random variables on the RHS of (3.19) are denoted by x, 

a: = log ^ ' I ' S 3.21 

{v\v) 

and the characteristic function of the identical distributions is 

F{k) = {e^'^PiO ■ (3.22) 



^If we replace the Gaussian distribution of C by a distribution where each element is real non-negative, 
uniformly drawn from the segment [0, 1], (3.16) does become an equality - in this case C has non-negative 
entries and the equality is a consequence of a Perron- Frobenius theorem; see chapter 8 in [24]. 
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The random variable on the LHS of (3.19) is denoted by y, 

{Vn\Vn) 



y = log 



{v\v) 



Its probabihty density is 



dk 
— ( 

27r 



-iky 



[F{k)T . 



(3.23) 



(3.24) 



We now expand x in e to order in the calculation of F{k) and assume that the 
expansion in e can be freely interchanged with various integrals. An expansion to order 
is assumed to be all that is needed, since an alternative treatment of the ensemble, based 
on a Fokker-Planck equation, would also need only expansions to order e^. 

By U{N) invariance we can rotate the vector v to point in the 1-direction, 



(3.25) 



Thus, the vector v has dropped out completely. We shall reuse the symbol v below. We 
now introduce some extra notation. 



Re Cii = u, Cji = Vj, Cij = Wj for j = 2, . . . , N . 
V and w are {N — l)-dimensional complex column vectors. We have to order 

X = 2eu + —{v^w + v'^w*) + e^w'^w . 
To calculate F{k) it suffices to know the distribution of u, v, w, 

P(n, w) = AA'e-^("'+"^'^+"'^"') . 
The integral giving F{k) is Gaussian and can be easily done, resulting in 

N-l 

m 



(3.26) 



(3.27) 



(3.28) 



e N 



1 ' 1/24 



(3.29) 



To the level of accuracy in at which we are working, we can write 



m 

The characteristic function of y is 



e iv e N 



-+m£2fc(l_ 1 ) 



Defining 



n 



(3.30) 

(3.31) 
(3.32) 
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the inverse Fourier transform gives 



P(^) = ^/^e-|^[*-^(^"^/^)^ (3.33) 



The Fiirstenberg theorems [15, 16] now tell us that almost surely 



1 , £2 / 1 



lim -log||T^„|| = - 1-- . (3.34) 

n^oo n Z \ JS J 

So far, and N have been kept fixed. We therefore conclude that for large enough n, 

||W^„||^et^'(i-i/^). (3.35) 

We now simply replace ne^ by the large number t and take N ^ oo, which is a relatively 
harmless limit. We conclude that 

7(t) ~ ^ for t ^ oo . (3.36) 

3.4.2 t^O 

We now wish to take e — > 0, keeping n and large but fixed. In terms of t, this would 
correspond to the asymptotic regime t — > 0. Since e is very small, we can try to expand to 
just linear order in e, keeping n and N finite albeit at large values. To linear order in e, 
the noncommutative aspect of the product is lost, and we can write 

Wn = ^i=i = e'^^ . (3.37) 

By an 0(n) rotation one can show that the matrix C = Sj=i Gaussianly dis- 

tributed. The distribution is fixed by its average and variance, 

{Cab) = 0, (ICabP) = ^ . (3.38) 

Here, 1 < a, 6 < are matrix indices. C is distributed exactly like in the Ginibre ensem- 
ble [23] . For iV — > oo we have 

spectrum(C') = {z; \z\ < 1} , (3.39) 

giving 

max{|u;|; w € spectrum (VF'n)} = = . (3.40) 

w 

We are therefore led to 

7(t) ~ for t ^ . (3.41) 
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3.4.3 -f{t) for all t 

Our subsequent work confirms the findings in [10] which, in turn, imply the existence of an 
annulus with a 'y{t) obeying our considerations. For = oo we shall find that the inverse 
function to 'j{t), which we call T^j) with 7 > 0, is given by 

r(7) = 27tanh|, T{-f{t))=t. (3.42) 

The previously presented asymptotic results are recovered. The two regimes, t ^ and 
t 00, differed in the order in e one goes to. With either choice, one obtains a finite 
expression if t is finite, so the truncation of the expansion in e is self-consistent. When the 
full limit n — > CO, e — > is studied at fixed arbitrary i, going to order should reproduce 
both asymptotic results in t, and we shall see that this indeed happens. Note that the 
crossover between the two asymptotic regimes occurs roughly where \/t = t/2, which 
means t = 4. It will turn out that as t increases the spectrum encircles the origin first at a 
critical value of t = 4. In some rough sense, this is the point where the lack of commutativity 
among the factors in the product of matrices becomes qualitatively important. It no longer 
is appropriate to think about the product as that of several matrices, each close to unity 
- a perturbative viewpoint. Some observables made out of the matrix product evolve in t 
smoothly, but others will develop a singularity. 

4. Average of products of characteristic polynomials 

For general N it is difficult to derive a closed formula for the distribution of all the matrix 
entries of W. We are interested in just the spectral properties of W. Even this is difficult 
to obtain for arbitrary finite N. Partial information about the distribution of eigenvalues 
can be obtained from the averages of characteristic polynomials related to W. These poly- 
nomials are generating functionals for various moments of the eigenvalue distribution. We 
shall denote the averages over the Cj hy (...) . The calculation of some simple characteristic 
polynomials is feasible. 

Obviously, (det(2; — W)) carries no information since expanding in z we see that only 
traces of products of M factors appear. The latter can be expanded in e and yield only 
factors of Cj , all of which vanish due to the invariance of P{C) under C Ce**. Hence, 

(det(z-t^„)) = (z-l)^. (4.1) 

The first non-trivial case is 

Q{z,z*) = {\det{z-Wn)\^) , (4.2) 

and from now on we focus on the calculation of the above in the limit. (By "the limit" we 
mean the limit n — > 00, e ^ with t = e^n held fixed.) 

If one applies large- factorization to (| det(z — Wn)P) (i.e., assuming that the average 
of the product can be replaced by the product of the averages) one gets holomorphic 
factorization, and all eigenvalues seem to have to be unity. For any f > 0, holomorphic 
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factorization will hold for z-values close to and z- values close to oo. These two regions 
are outside the annulus defined by 7(i). The full holomorphic factorized regime penetrates 
the annulus and will be connected for t small enough, but will split into two disconnected 
components for t larger than some critical value. There are two disconnected components 
when the eigenvalue support, always contained within the annulus defined by 7(t), fully 
encircles the origin z = 0. 



5. Saddle-point analysis of the basic random complex matrix model 

We wish to calculate Q as a function of t and see that at infinite N the transition we are 
looking for indeed occurs. As a first step we need some device to disentangle the nonabelian 
product defining W. Then, we can make the AT-dependence exphcit by integrating out all 
degrees of freedom whose multiplicity is iV-dependent. This allows us to take N ^ oo 
which, as usual, becomes a saddle-point problem. We analyze the saddle-point problem 
partially, only to the point where we can identify the transition we are after. 

5.1 Quark representation of characteristic polynomials for matrix products 

Let Xi,X2, . . . , Xn he n N X N square matrices of general structure. We are interested in 
the characteristic polynomial of 

W = XiX2---Xn. (5.1) 

We introduce nN pairs of Grassmann variables (quarks) : {'i/'j, V'j}j=i,...,n with the 
convention that V'n+i = i^i, etc. Let us define 

/n 
H [#,-#,-]e"' ^?=i ^^^^-^U . (5.2) 

i=i 

We claim that 

/„(Xi,X2, . . . .,X„) = det(y;" - W) . (5.3) 

We prove this by induction in n. For n = 1 the result is trivial. Assuming n > 2 we 
integrate over the pair V'nV'n to derive a recursion relation, 

W^'ln-l (^Xi, X2, . . ■ ., X„_2, = 4(^1. ^2, • • • , ^n) ■ (5.4) 

It is now easy to check that if the claim holds forn — 1 it holds also for n, and this concludes 
the proof. 

Obviously dct(ii;" — W) is invariant under cyclic permutations of the matrices Xj. In 
the quark representation this invariance is proved by "translating" the Grassmann pairs in 
their index j. 
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5.2 Making the dependence on TV explicit 

We start by defining a complex number a which depends on n such that 

2 = 6"-^, z = |z|e**, a = -\og\z\ + i- (5.5) 

n n 

with — vr < ^' < vr. We now introduce 4n Grassmann variables il)j,^j,Xj^Xj ^-^id write 

/n _ , _ 

(5.6) 

where ipn+i = and Xn+i = Xi- This can be rewritten as 

E"=i(e"V^,V'i+e-*XjX.)p-E"=iWiM,'/'.+i-X,+iArfx,) 



J] [di^jd^jjdxjdxj] e^^=i ^^■^^■) e 



i=i 

(5.7) 

We now perform several integration variable changes. We first switch the sign of Xj- 
Keeping the symbols x^X for t^is new variables we now replace Xj by Xj~i (with the 
understanding that xi ~^ Xn) and also replace Xj by Xj-i (again with the understanding 
that xi ~^ Xn)- When the integration measure for the new variables is written in canonical 
order, two (—1)^ signs cancel out. We are left with 



=E;=i(e"ViiV',+e"*X,X,)p-E"=i(V;iM,V'i+i+XiM+x,-l) 



\det{z-Wn)\^ = f Yl[di'jd'^jdxjdxj]e^^=''-'''^^^^'^'^'^''''^e- 

(5.8) 

The integral over the Ci , . . . , Cn now factorizes and can be done explicitly to sufficient 
accuracy in e to produce the correct t-dependent limit. The following equalities ought to 
be understood in the sense that they hold up to terms which vanish as n ^ oo, e ^ at 
t = e^n fixed. As usually is the case in stochastic calculations, we need to keep expressions 
correct to order e^, but not higher. We need to expand M only to linear order in e, as the 
next term in the expansion won't make a contribution after the C integration because of 
its phase invariance, and we end up with 

^^-{pMi,'-xM''x'^ ^ ^-4'i>'-xx'-§iix'xi>'-j^i^4''xx' _ (^^ g^ 

We used the external source formula (3.3) and included an extra minus sign obtained when 
a Grassmann variable was moved through an odd number of other Grassmann variables. 

We now separate the quartic Grassmann terms into bilinears by introducing scalar 
complex bosonic multipliers, Q and Xj, j = 1, . . . ,n, 

-N\(\^^-eiC^x'-Cx4>') ^ (5_10) 
-iV2|A|2g-e(AViV'-A*XX') (5.11) 

with 

iV2 

Mb = — . (5.12) 

TT 



e-#^x'x^' 




^-j^i'i^'xx' ^J^^ 


[ df^W 




N 


TT 
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The relative minus sign is needed to get the right sign in front of the quadrihnear Grassmann 
term. The integration measure is diJ,{() = dRe( dlm.(. 

The net effect was to replace the average over the complex matrix C by an average over 
the complex numbers (j and Xj, with all the "noise" now originating from the variables Q 
and Xj. This prepares the scene for making the dependence on N explicit. 

The limit we are after can be obtained from 



/n 
l[[di;jdi;jdxjdxjdii{Cj)d,,{Xj)] 



-wE"=ilC,P-^^E"=ilA,P 



X e 



X e 



- e E ?= 1 (Cj V-i Xj - 1 - C,* X J '/'j + 1 ) 



(5.13) 



One further change of Grassmann variables reduces the number of terms that are not 
diagonal in the index j, = tpj and Xj-i = Xj- Again, canonical ordering of integration 
measures leads to two (—1)^ factors which cancel out. Dropping the primes on the new 
variables we obtain 

/n 
]^[d^P,d^l;,dxJdxJd|,{CJ)d^i{X,)]e-''^U\^^\^^^^ 

3=1 



X e 



Carrying out the integral over the Grassmann variables we get 

{\det{z-Wn)\^)=M2K j i\W{Qdti{Xj)]e-''^UM'-'''^^^^^ 
where 



N 



y4 = e^r^-l-eA, D = e^*T-l + £A^ B = -eZ , C = eZ'^ 



(5.14) 



'A B\ 

cb)^ 

(5.15) 
(5.16) 



with 



/O 1 
1 


\1 



o\ 



1 
0/ 



Z = diag(Ci, • • • , Cn) , A = diag(Ai, . . . , An) • 



(5.17) 



Using known formulas on determinants of block matrices, we have 

det j^^J =det(AL'-^C^-^5) =detAdetDdet(l + £2zU-iZD-i) . (5.I8) 
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5.3 The trivial large- AT saddle and its domain of local stability 

Since the A^-dependence of the A-integral is of the form 



kJ n[dMA.)]e-^^^"-l^^l^det^ ^] , (5.19) 



we evidently get a trivial saddle point Aj = 0, j = 1, . . . , n for large N due to the dominance 
of the iV^ term. In this limit, we can therefore focus on the remaining (^-integration with 
the replacements 

Ao = e''T^ -1, L> ^ Do = e'"*r- 1 = ylj , (5.20) 

which yields 

{\detiz-W^r)^M:l n[rfMO)]e-^^"-l^^l^det^ (^^° j^^ . (5.21) 

This expression is exactly equal to the result which we would obtain without restricting 
the determinant of Wn- Integrals over complex A- variables were needed to decouple quartic 
Grassmann terms arising from the second term in the exponent of (3.3), which simply does 
not occur in (3.7). The boundary of non- vanishing eigenvalue density is therefore equal for 
Wn G GL{NX) and Wn G SL{N,C) in the large-iV limit. 

Using again the notations A and D for Aq and Dq, we have 

det ^ = det{AD - ACA-^B) = \ det A\'^ det(l + A^^ Z{A^)'^) . (5.22) 
The matrix T implements cyclical one step shifts and obeys T" = 1. Hence we can write 

n 

A-i = — V = ye-^"^^ (5.23) 

s=l 

Each entry in A~^ gets exactly one contribution from a single term in the sum over s above. 
A~^ is a circulant matrix, which means that it has identical entries on lines parallel to the 
main diagonal and periodic with period n on lines parallel to the anti-diagonal. 

The large- limit will obviously lead to saddle-point equations which will be satisfied 
at Q = since the ^j, Q enter only bilinearly in the integrand. If this saddle dominated at 
infinite A^, we could replace Wn by a unit matrix, 

{\det{z-Wn)\') = \z-l\^^ . (5.24) 

This means there are no eigenvalues at any z ^ 1 in the complex plane. In particular, the 
eigenvalue density in the complex plane, scaled to be finite at infinite N, is zero everywhere 
(except at the singularity z = 1). We shall refer to this saddle as saddle A. Where saddle 
A gives the correct answer, (I det(z - Wn)\^) = \z- Ip^ is the absolute value square of a 
holomorphic function in z and there is no finite eigenvalue surface density. An eigenvalue 
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surface density will develop in regions of the complex plane where saddle A is displaced 
by another saddle, saddle B, which destroys holomorphic factorization. To determine 
where saddle B must take over we find the regime where saddle A is no longer locally 
stable. Comparison with numerical simulation shows that saddle A is always dominating 
whenever it is locally stable and that saddle B indeed produces non-zero surface charge 
density. Thus, at the boundary of the domain of stability of saddle A one has a continuous 
transition to regions with non-zero surface eigenvalue density. We do not calculate saddle 
B explicitly. 

5.3.1 Determination of the boundziry of the domain of stability of the trivial 
saddle 

To determine the domain of local stability of the trivial saddle point we need to calculate 
the Gaussian form of the integrand around saddle A. To quadratic order in Cj^Cj we have 



detM^ =|det(l-e'^Tt)|V^^. 



(5.25) 



We are interested only in F, 



F = tr ZU-'Z{A^)-' = ^c;|(^-i)^,|2^, = E^l^O ■ 



(5.26) 



Indices are understood modulo n, with the fundamental interval taken from 1 to n. The 
matrix K is also circulant, so its eigenvalues are determined by the discrete Fourier trans- 
form of the last row of the matrix K which defines the entire matrix in terms of an n-term 
series, Knj = Kj, with 



Kj = 



1 



1 



-j{a+a*) 



(5.27) 



The eigenvalues of K are 

1 2 n 



1 



1 — e " "e 



fep-((T+0-*) 



/_ 1 _*\ 27vik 



(5.28) 



Going back to the original variables, the condition for local stability is 



IP 



Re 



2 2-Ri 1 
\z\ne n 



< 1 



(5.29) 



for all = 1, . . . , n. It is easy to see that if > 1 the inequality is strongest for the k = n 
case; the same is true for \z\ < 1. Hence, the condition holds also for all A; < n if it holds 
for k = n. We end up with a determination of the region of local stability of saddle A, 



Taking n ^ oo at constant t in Eq. (5.30) gives the chargeless region, 



f IzP — 1 

1 > 91 112 'i I 1 ^ (5-31) 
2|z — log \ z\ 



in agreement with Eq. (83) of [10]. 

It is easy to see that the points on the boundary, separating charged and chargeless 
regions, having maximal or minimal absolute values are located on the positive real axis. 
This means that the function 7(t), defined in Sec. 3.4, has to fulfill 

, , t eT(*) + 1 , , 

which is equivalent to Eq. (3.42). 

Note that the exact invariance under inversion and complex conjugation of z has been 
restored in the limit, although it was lost at finite n because of the truncation in the 
expansion in e to second order (which was all that is needed to get the correct limit). 
Therefore, as explained earlier, one can look for the transition point by just focusing on 
the unit circle. The portion of the unit circle which resides in the chargeless region is 

^<|z-l|^ 1^1 = 1. (5.33) 

For t < 4 there is an arc centered at z = —1 which resides in the chargeless region. The 
end points of this arc are at the the two angles satisfying cos(^) = 1 — t/2. When t > 4 
the charged region contains the unit circle and hence becomes multiply connected. The 
last point to be engulfed by the charged region as t increases is the point z = —1. 

5.3.2 More detailed study of the neighborhood of the critical point 

To better focus on the shape of the boundary on both sides of the transition point t = 4, 
it is useful to employ the following maps, 

z = and z = oo map into u = ^1/2, and z = 1 maps into u = oo. There is always 
charge at z = 1, so in the n-plane the charged region extends to infinity. The \z\ = 1 
circle maps into the imaginary axis in the u-plane, and z = —1 maps into the origin u = 0. 
Inversion in z becomes u — > —u. The real- 2; axis maps into the lesl-u axis. The region 
{Imz > 0} n {\z\ > 1} maps into the region {Imn < 0} H {Reu > 0}. Reflection about 
the real axis {z z*) maps into reflection about the real axis in the u-plane {u u*), 
and reflection with respect to the unit circle {z 1/-2*) corresponds to reflection about 
the imaginary axis in the n-plane (n — > —u*). Our problem has these symmetries in the 
z-variable, so it suffices to analyze one of the four quadrants in the n-plane and get the 
result for the other quadrants by reflection through common axes in the n-plane. 
We have seen that eigenvalues are restricted to the annulus 

e-^W < \z\ < e^W . (5.35) 
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Therefore, the complement of this annulus is contained in the chargeless region. Under the 
map, the two circles \z\ = e^'^'-*^ go into two circles with non-overlapping interiors in the 
tt-plane. 



^coth(±7(t)) 



1 



(5.36) 



(2sinh7(t))2 ' 

The eigenvalues are restricted to the exterior of these two circles (the image of the annulus) 
and therefore the chargeless region (5.31) includes the interior of these two circles. 
Denoting Re u = Ur and Im u = Ui, the chargeless region is found to be 

0<uf < Ur coth(tur) - - 1/4 . (5.37) 

For large enough the above inequalities self-contradict, showing that the chargeless 
region is bounded in the n-plane. 

The chargeless region in the positive quadrant of the u-plane determines the entire 
chargeless region by reflections from quadrant to quadrant through a common axis. When 
Ur ^ we have 

< ^ - ^ + - l) + 0{ut) . (5.38) 

For i < 4 there is a portion of the imaginary n-axis inside the chargeless region; hence, the 
charged portion of the imaginary n-axis has a break around the origin, and this maps into 
the unit circle in the z-plane having a gap around z = —1. For t < 4 there is one connected 
chargeless region containing the two circles from Eq. (5.36). For t > 4, Ur = is not in the 
chargeless region. Consequently, the entire imaginary axis is in the charged region. The 
chargeless region is split into two disconnected regions, each containing exactly one of the 
two circles of Eq. (5.36). 

Exactly at the transition, when i = 4, the boundary in the vicinity of the origin is 
given by the two lines 

Ui = ±—^Ur + . . . (5.39) 
v3 

This critical contour is a slightly deformed figure-8, laying horizontally along the leal-u 
axis and symmetrically with respect to the imaginary-u axis. The midpoint of the 8, which 
resides at the origin of the u-plane, separates along the real u-axis as t is decreased from 
t = 4, splitting into two disconnected regions, and separates along the imaginary-u axis as 
t is increased from t = 4, becoming one single connected region. 

5.3.3 Connection to the inviscid Burgers equation 

The formula (5.31) for the boundary of the chargeless region leads us to introduce the 
following map from the complex plane onto itself, 

Z(u,t) = ^|e-*" . (5.40) 



u 2 



For Reu ^ Eq. (5.31) is equivalent to 



, , . , f > 1 for Re u(z) > , , , 

\Z{u,t)\ { V ^ ' 

I < 1 for Re u{z) < , 
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and the boundary between the chargeless and charged regions is given by 



\Z{u,t)\ = l, Reu/0. (5.42) 

For Reu = 0, the boundary is found using (5.38). 

Only at t = is Z{u, t) one-to-one. For nonzero t, Z{u, t) has an essential singularity 
at ti = oo which prevents an analytic definition of an inverse, U{z,t). One is therefore led 
to look for a local definition of the map, by a partial differential equation. We differentiate 
the equation 

Z{U{w,t),t) = w (5.43) 

with respect to w at fixed t and with respect to t at fixed w. We find then that U{w,t) 
obeys 

= ^"9.; • ^'-''^ 

This is the inviscid complex Burgers equation (the Hopf equation), up to a change of 
variables w = e~^. This equation plays a central role in two-dimensional YM and in 
the equivalent multiplicative random unitary matrix model [3]. It comes with the initial 
condition 

U{w,0) = U'^)^u{w) = -lcoth^. (5.45) 



Defining ^{w,t) by 



Ui-^,tHl/2 f/(.,t) = lcothii^, (5.46) 

U{w,t) - 1/2 ' ^ ' ^ 2 2 ' ^ ^ 

we find that (5.43) is equivalent to 

e«(-'*)-i™*^^=«;. (5.47) 

This can be viewed as an equation for ^, which in turn defines the solution of the partial 
differential equation with the desired boundary condition. The variables x and ^ are, up 
to some factors of zizi, identical to the variables denoted the same way in [3]. 

Equation (5.42) identifies the boundary separating the charged region from the charge- 
less one with the image of the circle {wl = 1 in the u-plane under the map u = U{w,t) for 
Re n 7^ 0. It is well known that as t increases from zero, depending on the initial condition, 
singularities can be generated at finite t > 0. In our case, we have seen explicitly that at 
t = 4 a singularity is generated. 

5.4 Precise relation to the model of Gudowska-Nowak et al. 

Instead of (3.4), we could define 

Mj = l + eCj . (5.48) 

This is the exact form of the model of [10] . It looses the inversion symmetry in z at finite n. 
With this definition our formulas become exact even for finite n. The limit n — > oo, e — > oo 
with t = ne^ held fixed will not change if (3.5) is replaced by (5.48) and the inversion 
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symmetry is recovered. However, with (5.48) we can explicitly work out a few low-n cases 
and test them either numerically or by more direct analytical means. We did this only for 
n = 1,2. For these values of n we obtained agreement with [10], providing an additional 
check on our method, which relies only on information captured by the observable Q{z, z*). 

This "linear" model (as opposed to the previous, "exponential" version) is more con- 
venient for numerical simulations, as one does not need to exponentiate a large matrix. 

5.5 Numerical results 

In general, we have not determined the global stability of our saddle, nor have we identified 
explicitly the competing nontrivial saddles. Therefore, we need a bit more evidence to 
establish the transition. We do this numerically. Numerical checks were first carried out 
in [10], and we confirm their results. 

We want to work out the inverse map U{w, t) explicitly. First, we agree to focus on one 
quadrant and use the symmetries to get the solutions in the other quadrants. We choose 
the quadrant Ur,Ui > 0, where short hand for the real and imaginary parts of 

the function U{w, t). Similarly, we denote the real and imaginary parts of v = 5^)^ by Vr 
and Vi respectively. We are given Vr,i and need to find Ur^i- The equation to be solved is 

"-±lIle-'- = ^='-±^^ (5.49) 

u-1/2 v-1/2 ^ ' 

Here t is real but can have any sign, u = Ur + iui and v = Vr + ivi. Taking the absolute 
value we get 

|np + 1/4 = n.^coth(6 + tti^) , (5.50) 

where 

tanh6 = , ,0 — rr • (5.51) 

|t'P + l/4 

This gives Ui in terms of Ur- The boundary of the chargeless region corresponds to Vr = 0, 
see Eq. (5.37). 

A second relation is obtained by looking at the phase. We first write 



1/2 



V - 1/2 

Then 



v + l/2 



1/2 



n + 1/2 
u-1/2 



n+ 1/2 



u-1/2 



(5.52) 



a = f3 + tui , (5.53) 



where, by convention, |q!|,|/3| < vr. /3 is known, and we have an expression for 
function of Ur and other known quantities from above, a is a function of Ur^i so the above 
is a transcendental equation for Ur (or a), 

I up — 1/4 Ui Ui 
cos a = T-R — TTT , sm a = — r-^ -— - , tan a — 



|u2-l/4|' |u2-l/4|' |u|2-l/4' 
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The simplest form of the transcendental equation is obtained using the formulas for the 
tangent, but the other equations are needed to resolve some discrete ambiguities. 

All numerical simulations were performed with n = 2000 and = 2000 for ensembles 
consisting of about 500 matrices (without a restriction on the determinant). In the figures 
below, the solid lines correspond to the analytically derived boundaries, which are in very 
good agreement with numerical data. Figure 1 shows eigenvalue distributions in the z- and 
li-plane for t = 3, t = 4, and t = 5. Numerical tests confirm that the topological transition 
of the domain of no n- vanishing eigenvalue density occurs at t = 4, when the domain 
becomes connected at z = —1. This corresponds to the imaginary n-axis completely lying 
in the domain of eigenvalues. Note also that the eigenvalue density in the u-plane is indeed 
symmetric under reflections at the real and imaginary axis, which is related to inversion 
symmetry in z. 

Figure 2 shows eigenvalues for t = 12 and affirms that for large t the boundary ap- 
proximately consists of two circles with center at z = and radii exp(ibt/2). 



6. A generalized Gaussian random complex matrix model 
6.1 Definition and general properties 

The basic complex matrix model can be generalized to interpolate between the case in 
which the individual factors are unitary and the case in which they are hermitian. Writing 
the matrix C of each factor as C = Hi + iH2 with Hi^2 hermitian and traceless, we take 
the probability distribution of Hi^2 to be 

P{Hi,H2) =Afe~^^^''"'^^''^'^ (6.1) 



with u>i^2 > 0. This produces the following correlation functions 



(CijCik) = {C*jCik) = ^(wi - uj2)5ik5ji - -^{^1 - ^2)kj5ik (6.2) 



4,1, , . . 1 



and 

{CijCjk) = ^('^i + ^2)SikSji - -^{^i + ^2)kjkk ■ (6.3) 
As a consequence, we have 



[e^^e^^X^ = (l + 2-1 (l -^)e' + 0(e^)) 5,- . (6.4) 

Using the same analysis as in the original model we learn (see below) that the spectrum 
of Wni in the limit n — > oo, e — > with t = e^n held fixed, will now be restricted at any 
to the annulus^ 



*Now it would be meaningful to make e complex, since for uji ^ UJ2 the phase invariance C e'*C has 
been eliminated. One still has invariance under a sign switch of C and only enters. A complex e would 
produce a complex t. 



- 24 - 



6 
4 
2 

I " 
-2 

-4 

-6 

10 




-1 1 



4 5 6 7 8 



-10 

15 
10 
5 

I " 

-5 
-10 
-15 













\ 

\ 




■ j 

/ 











2 4 6 8 10 12 

Refz) 




10 15 20 



Refz) 



0.2 



-0.2 



-0.4 









X 


































1 


-0.5 





0.5 








Re(ul 







-0.5 





Rc(u) 



a 




0.5 




Figure 1: Scatter plot of the eigenvalues of W in the z-plane (left) and in the u-plane (right) for 
t = 3 (top), i = 4 (center), and t — 5 (bottom). 



where the function 7 is defined in Sec. 3.4 and fulfills Eq. (3.42). As wi — > the spectrum 
is restricted to the unit circle \z\ = 1, and we are in the unitary case. For = a;2 = 1/2 
we are in the original model. For a;2 = we are in the hermitian case, which produces a 
complex matrix upon multiplication of the individual factors. 



6.2 Large- AT factorized average 

For uJi ^ u)2, 

(det(z - Wn)) = J{z) 



(6.6) 



- 25 - 




Figure 2: Plots of the eigenvalues of W in the z-plane for t — 12. The plot on the right is an 
enlarged version of the plot on the left. 

is no longer equal to {z — 1)'^ , but to a more complicated polynomial in z. The polynomial 
J{z) is completely determined by the two-point function of the matrix C. 
Where factorization holds we have 

{\det{z-Wn)\') = \J{z)\' (6.7) 

and no finite surface charge density. 

The polynomial J{z) can be read off from previous work [2] on a product of random 
unitary matrices. Take uJi < uj2- The correlation functions in Eq. (6.2) tell us that J{z) 
depends only on the difference uj2—uji, so we could simply set wi = 0, and then we obviously 
are in the unitary case with C = iH2. So, for uj2 — loi > Q the dependence on L02 — can 
be absorbed in t by a rescaling, 

t± = t(u;2 iwi), t+>|i_|. (6.8) 

From [2], we get for the SU{N) case 

\ — / A'^ \ AT I 1 t'k(N-k) 

J{z) = QN{z,t^) = E U ^ (-1) e"^^ ' (6-9) 

fc=0 ^ ^ 

where t' = t_(l + 1/A^). For large N nothing is lost by ignoring the difference between t' 
and t-. From now on we set = t_, 

N 

Sr-^ ( N\ AT I 1 t_k{N-k) 

•^(^) = E J z^-\-l)''e ^ . (6.10) 

k=0 ^ ^ 

We also have an integral representation, 

j(,) = \l^ r dX e-f \z - e-*-(^+V2)l ^ . (6.11) 
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In [2] it was shown that for i_ > the above polynomial has all its roots on the unit 
circle. This was done by interpreting the integral as the partition function of a classical 
ferromagnetic spin 1/2 model in an external magnetic field determined by z. 

For a;2 — < we need to analytically continue to negative t_ . As observed in [2] this 
is evidently possible in the polynomial form. However, it no longer is true that all zeros 
are on the |z| = 1 circle. One can also analytically continue the integral expression, 



j(^) = r dAef*-^' L_e-*- (^^+1/2)1^, (6.12) 
V 27r y„oo L 

where, by Eq. (6.8), now t_ < 0. 

Equations (6.11), (6.12) can be combined into one line-integral expression, 

j(,) = / d\ e-f \z - e-*-(^+V2)l ^ . (6.13) 



27r 



Here, t_ is real of either sign and £ is the real axis (from -co to oo) for t_ > and the 
imaginary axis (from —ioo to -|-ioo) for t_ < 0. For t_ > we take \ftZ > 0, and for 
t_ < we take ^JtZ = —i^—t^ with \/—t^ > 0. 

6.3 An exact representation of the generalized Gaussian model 

To determine for which values of z the factorized formula no longer holds and one expects 
non-zero surface charge density as a result of the loss of holomorphic factorization we 
need the analogue of Eq. (5.21). The method of getting at this formula is the same; the 
one complication is that in addition to the complex noise factors ("j and Aj one needs to 
introduce two additional real noise factors and dj. These extra noise factors are needed 
because more quadrilinear Grassmann interaction terms need to be decoupled. Still, at the 
end the dependence on N is made explicit. 

It is also convenient to introduce the notation 

(j+ = ^e\ujx^^2) , (6.14) 

f v^|e2(u;2 - u;i)| for W2 > wi , 

LU- = < (6.15) 

[i^\e'^{LL>2 - uji)\ for tJ2 < '^i • 

Since the integral over the variables Xj is again of the form (5.19) it can be trivially 
approximated by a saddle point at the origin in the large- limit. We are then left with 

/n 
n [dfiiCMjdOj] e-^^^=i(l^^l'+^«'+i^')det^ 
1 

where A/'c = N/2tt and 




^ = - ( 1 - ( 1 - ^ ) - ^_ y^l + 1h 1 + e^Tt, (6.17) 



= - 1 - ^cuM 1 - 4- ) - u;_ Jl + i-e + e'^Tt, (6.18) 



2 " V V A 

B = -Lo+Z = -C7t (6.19) 
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with 



/O 1 
1 


\1 



o\ 


1 

oy 



(6.20) 



and 



Z = diag(Ci, . . . ,Cn) , S = diag(^i, . . . , G = diag((9i, . . 
6.4 Region of stability of the factorized saddle 

The identity 

= det{A) det{D) det(l - D^^CA-^B) 



9n)- 



(6.21) 




(6.22) 



shows that the variables C only enter as bilinears CjCk- large N the ^-integral will be 
dominated by some saddle point, and Cj = is a trivial solution to the C saddle-point 
equations, for any Q, H. Where this saddle point dominates, z is in the chargeless region. 
The reason is that we have C = i? = at this saddle, and then the remaining H and Q 
integrals factorize and we have 



det(z - Wn 



n [d^jdd,] e-^S"-(i«M^.')det^(ylD) 



/n 



(6.23) 



Since 



det(^) = (-1 



n 1 

j=i \ 



-UJ_ 



1 

iV2 



1 



1 + ^e. 



(6.24) 



depends only on z, and not on z*, we have holomorphic factorization. Therefore, it is 
just the structure of the (^-saddle which determines that z is in a chargeless region. The 
holomorphic factor in Eq. (6.23), however, is needed to determine local stability. At the 
C = saddle, from [2] we know that in the e — > limit we have 



/n 
n [^^^-l e-''^"-^('^^")det''{A) = (-l)^("-i) J(z) . (6.25) 



In the large- limit the simplified formulas (6.11), (6.12) apply. 

To determine local stability we need the saddle point which dominates the integral 
over the in (6.25). The derivation leading from (6.25) to (6.11) shows that for the case 
of tLi2 > uJi, at the saddle one has = ^, where 



A. 



(6.26) 
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with finite A^. To leading order in N we then have 



det(^) ^ (-1) 



n-l 



j=i \ 



2n 



{-ir-'[z-e 



-t_(A,+l/2)i 



(6.27) 



As is a saddle point of the integrand of (6.11). With obvious changes, a similar story holds 
for UJ2 < wi. 

We end up with the following expression for the matrices A and D needed for the 
analysis of the stability of the trivial saddle under variations of C, 



A 



l--(As + l/2) 
n 



1 + e'^T'l' 



Here A^ has to satisfy 



A,, 



1 



(6.28) 



(6.29) 



^ 2gt_(A,+l/2) _ I 

The appropriate contour C in Eq. (6.13), whose endpoints at infinity are fixed, will be 
deformed to A^, and we assume that the integral will be dominated by one single saddle 
point as long as one is in the charge less region. 

It is now convenient to define ii = Xg + 1/2. Then the map from w to z acquires a 
simpler form. 



z = z{u) 

In terms of the map U{w, t) we have 



u + 1/2 
■u- 1/2' 



-t-U 



Z{u, t- 



u = U{z,t_) 



(6.30) 



(6.31) 



where we also allow for t_ < 0, corresponding to "backward evolution". 
The stability of the trivial saddle is now determined by u from 



det 



1+ujIz^A~^Z{A^)-^ 



with 



A 



1 + e'^rt = e" 



(6.32) 



(6.33) 



where, as before, e"'^ = z. We can drop the prefactor e " because in the determinant 
there is an extra prefactor. We therefore end up with 



det 



with 



and 



1 + *±zti-iz(it)-i 

n 

A = -l + e^T^^ 



u- 1/2 



(6.34) 

(6.35) 
(6.36) 
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z is defined by the complex number ii which should be determined by z and t_ in (6.30). Let 
the unit circle = 1 be parametrized by |s| < vr, with u; = e**. For \z\ 7^ 1, the boundary 
separating the charged and chargeless region is defined in the z-plane by z = /(s), given 



Note that > For t+ < 4, the boundary intersects the unit circle in the z-plane at 
the points 



In complete analogy to the basic model, setting aside the restriction det(Wn) = 1 does not 
affect the boundary in the infinite- limit. 

6.5 Numerical results 

As mentioned above, the linear model is much more convenient for numerical simulations 
than the exponential one. Performing a similar stability analysis for the linear model, it 
turns out that the boundaries of the domains with non- vanishing eigenvalue density for the 
two models are equivalent up to a scaling by a factor of exp(— 1_/2). The following figures 
show perfect agreement between numerically obtained eigenvalue domains and analytically 
determined boundaries for the linear model (data points as well as predicted boundaries are 
scaled by the corresponding factor of exp(— 1_/2)). Therefore, we expect that the stability 
analysis gives the correct boundary for the exponential model, too. 

Figure 3 shows results of numerical simulations for wi = 1/10, u;2 = 1/2, with all other 
parameters as in Sec. 5.5. The topological transition occurs at t = 20/3, which corresponds 
tot-)_ = ^(2+]^) = 4 in agreement with the prediction. 

Figure 4 is generated with uj2 = 1 and uJi = 1/2000 for t = 1. Since this is already 
close to the unitary model, the eigenvalues are restricted to the vicinity of the unit circle 
in the z-plane, which corresponds to the imaginary axis in the u-plane. As t+ is below the 
critical value we get no charge around 2; = — lorn = 0. 

7. Beyond infinite N and the associated saddle-point approximation 

As mentioned in the introduction, we re-analyzed the model of [10] because we guess that 
it is a universal representative of the large- phase structure of large classes of complex 
matrix Wilson loops. We shall refer to this hypothetical class as the large- universality 
class. 

Like in the unitary case, we would like to study in more detail the approach to infinite 
N and see what the matrix model universal features of this transition are. To do this 
we need a more convenient finite- representation of the average of the product of two 



by 



/(s) = Z(C/(e*^^+),^_) for Re C/(e*^ t+) / . 



(6.38) 




(6.39) 
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Figure 3: Scatter plot of the eigenvalues of W in the z-plane (left) and in the u-plane (right) for 
LUi = 1/10, ^2 = 1/2 and t = 5 (top), i = 20/3 (center), t = 10 (bottom). 



characteristic polynomials we have been looking at. More precisely, we would like to first 
take the continuum limit £ ^ without making any assumptions about how large N is, 
and only later take N large. We turn to an outline of how one could achieve this; however 
we have not carried this to completion because, as will become clear, a full analysis keeping 
full exact A^-dependence is complicated. We hope that using our technique one could learn 
which subleading terms in the large- limit can be dropped without changing the universal 
properties of the approach to the large- limit. We start with the unitary case where the 
problem has been solved in [2] and present an alternative way of deriving that solution. 
This alternative way has the potential to generalize to the complex matrix case. We first 
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Figure 4: Scatter plot of the eigenvalues of W in the z-plane (left) and in the u-plane (right) for 
uji = 1/2000, UJ2 = 1 and t = 1. 

look at the simplest complex matrix model and then at the more general one. 
7.1 The unitary case 

As before we consider the SU{N) case but drop some irrelevant 0{1/N) corrections to 
keep the formulas simple. When considering products of unitary matrices it suffices to 
look at the average of the characteristic polynomial, and there is no need to calculate 
the average of its absolute value squared. This is a significant simplification. Employing 
a quark representation (for details see [2]) we found the following representation of the 
average characteristic polynomial. 



(det(z - W)) 



n 

i=l 



NdX 



N 



11(1 -y-^^^-: 



e 2 



AT v^n \2 



(7.1) 



Here, W is the product of n unitary matrices all close to the unit matrix. Introduce the 
notation 



p{X)dX 



N N 




e^^^ dX 



and 



n 



eXj) 



(7.2) 



(7.3) 



We are only interested in the limit n — > cx), e — > with t = ne^ kept fixed. We eliminate 
the A-independent term without affecting the limit by introducing new variables 



e 



(7.4) 



We now proceed by finding the probability density distribution for the variable cin- 
In other words, we look for a way to perform the integral over all A- variables keeping the 
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product we are interested in fixed at some arbitrary value a. This can be done in the hmit 
we are interested in. There, the probabihty density for a would be P{a;t)da. P{a;t)da 
is obtained from the Pn{d)da, the probability densities governing the variables a„ at step 
n. The basic step is to derive a recursion relation for the Pn{o) and take the limit on 
that recursion relation. In the limit the recursion relation becomes a partial differential 
equation of first order in t and second order in a. The reason for being first order in t is that 
Pn{a) is determined by P„_i(a); it does not depend on Pk{a) with k < n — 1. This reflects 
a "Markov property" of the multiplicative structure of the original matrix ensemble, and 
hence is a fundamental property [14] . The partial differential equation is of Fokker-Planck 
type. In conjunction with the boundary condition lim(^o+ P{^] t) = ^{o, — 1) the equation 
determines P{a; t) completely. 
Suppose we are at step n, 

a„ = o„_i(l - eXn) . (7.5) 

We now drop the subscript n, i.e., = a, A„, = A and a„_i = o + 6d. To the order to 
which we need to keep terms we have 

da = -a+ = a{eX + e^X^) . (7.6) 
1 — £A 

Probability logic tells us that the probability to find a„ between a — ^ and a + ^ is given 

by 

Pioh(^a-^ < dn < a+^^ = J dAp(A) Prob(^a-^ < a„_i(l-eA) < a+^^ . (7.7) 
Probability densities depend on the measure, 

Probfa - ^ < a„ < a+ ^) = P„(a)Aa. (7.8) 



2 2 

We end up with the recursion 



Pn{d) = j d\p{\)Pn^i{d + 5d){l + e\ + e^\^) , (7.9) 

where the last factor comes from the change in measure. We now expand to second order 
in e and do the A-integral, 



-2 



Pn{d) = Pn-i(a) + ^ 
We finally can take the limit. 



„ . _aP„„i(a) , l^^d'Pn-i{d) 



Pn~i{d) + 2d — h -d 



dd 2 da^ 



(110) 



Ar^=P(a;*) + 2a^ + la=?^. (7.11) 

ot oa 2 aa^ 



This equation can also be written in the form 

1 



NdtP = daidP) + ^-dad^daP , (7.12) 
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showing explicitly that the integral J Pda is time-independent and therefore equal to unity, 
its initial value. 

The equation with the delta function initial condition has the following solution [14], 
describing a log-normal distribution. 



We can now write an expression for the average characteristic polynomial, 

(det(z-Ty))= [ da(z-ae-^Y P{a;t) . (7.14) 



We change integration variables, 

a = e~*''-2^, (7.15) 



to finally obtain 

{det{z-W)) = / dfilz-e^^-^f"-^)" e'"^ . (7.16) 




After dropping the t/2N term in the exponent inside the parenthesis in the integrand we 
obtain the result of (6.11). 

One does not need to solve the Fokker-Planck equation exactly in order to get the 
large- limit, because N plays the role of 1/^ (with Euclidean time though) in (7.12). 
At large A^ the Fokker-Planck equation reduces to a Hamilton-Jacobi equation for S with 



2 



N— = -a^ — . (7.17) 
dt 2 \dd J ^ ' 

The solution that satisfies the initial condition is particularly simple, 

5 = -^log2a. (7.18) 

As usual, there is a prefactor to which contains additional t-dependence, but this factor 
is not needed for the large- A limit. There are a few other terms which can be ignored in 
the large- A limit. 

In the expression for the average characteristic polynomial we now have two terms that 
are exponential in A^, 

The saddle-point equation is 

_ j. 

e 2 lloga 

^ (7.20) 



z — e 2d ^ ^ 

-t\ 



Define now d = e and rearrange the saddle-point equation slightly to obtain 

A = (7.21) 

1 - ze*^+2 
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This reproduces the saddle-point equations we had before. 

The main conclusion is that, as it often is the case in the context of large- models, 
one has a "quantum" -like equation for finite A^, with 1/N playing a role analogous to h. 
The large- limit is then "semiclassical" , with the "quantum" equation being replaced by 
a classical one, in a variant of the WKB method. 

This is what we would like to duplicate in the complex matrix case. 

7.2 The basic product of random complex matrices 

7.2.1 An exact map to a product of random 2x2 matrices 

In the following, we do not restrict the trace of Wn to keep the analysis as simple as 
possible. As explained in Sec. 5, we can then start from Eq. (5.13) without the A-integrals, 

/n 
n[dV^,#,^ix,rfx,^^^(0)]e-'^^-^^^''^^-^^-^-^^*^^'^^+^^ 

X e" E"=iW,V'i+i+x,x,-i)gE;=i(e"V;,V',+e'^*x,x,)g-^E"=i IGf _ (7_22) 
We now change notation by introducing two-component Grassmann variables ^>j, 

= (^^^^^^ , = (Vi, Xj) ■ (7.23) 

Equation (7.22) can then be written as 

„ n 

(I det(z - Wn)\^) = (-l)^("-i) J J][d$,-d$,-d//(0)]e-^^^=^ '^^l' 



X exp 



.i=i 



-I I ' ' \eq 



(7.24) 



where the factor (— l)^("~i) results from bringing the integration measures in canonical 
order after the index shift in Xy 

(e'^ —eC,\ 
^ ^ \. Now, using the 

identity (5.3), we can integrate over all the Grassmann variables, obtaining an expression 
containing another random matrix product, but this time the matrices are just 2x2. 

(|det(z-l^„)|2) = (-z)^|n[dM0)]e-^^^-^l^^^l^^^^ l^l-^n^jj ' (7-25) 



where 
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Absorbing some phase factors into the variables Q , the matrices Qj can be further simpHfied 
to 



V eq e 2 ; 

The matrix after the prefactor is an SL{2,C) matrix. Hence, the product of Qj matrices 
will be, up to a multiplicative factor, also in SL{2,C). The new multiplicative random 
matrix model defines a stochastic process on the SL{2,C) manifold. Unlike the N x N 
matrix ensemble, the 2x2 matrix ensemble has no inversion symmetry and only a restricted 
conjugation symmetry. 
Finally, we write 

/n 
n[^/^(0)]e-^^"- l^^l'det^ (1 - e-*A„) . (7.28) 

Here, A„ = HLi and 



,^,^,N-V.a..= |,n^.0^j (7.2.) 



The new multiplicative matrix ensemble has detlj = 1, so A„ is restricted to SL{2,C), 
which in turn implies 

/n 
ll[dfx{Cj)]e-''^"=^ l^^l' (tr A„ - 2cos*)^ . (7.30) 

i=i 



In the new ensemble we have invariance under complex conjugation and conjugation by a 
U{1) subgroup, 

e-^^> J M e'^J 



n .-ioAYA . D- (7-31) 



Therefore, the probability density of A„ will be invariant under 



So, the FP equation for this case, 

^ = @N^N , (7.33) 

is an equation for a function Sjv of six real variables including t {z is a parameter). Qn 
is a linear partial differential operator of second degree in five real variables. There are no 
terms from the measure if we pick the latter to be 5*^(2, C) invariant. 

The dependence on N is explicit. All second-order derivative terms carry a 1/N 
factor. Among the first-order derivative terms some have a 1/N factor and others are 
A'^-independent . 
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One can write down an exact integral expression for Q{z, z*) in terms of Sjv(A„; \z\, t), 

Qiz,z*;t) = \z\^ j dn{An)^NiAn;\z\;t){tiAn-2cos^f . (7.34) 

In the infinite- limit, Q{z, z*) would be given by a dominating saddle point, and sublead- 
ing corrections would identify the relevant large- universality class. 

An analysis of the Fokker-Planck equation for arbitrary z seems too complicated to 
attack directly, so we restrict our attention to the unit circle, \z\ = 1. 

7.2.2 The case of |2| = 1: Simplifications for large N 

We now set \z\ = 1. This simplifies the operator Qn, eliminating the first-order derivative 
terms that had no A^-dependence. Consequentially, at \z\ = 1 one has 

= , (7.35) 

where G has no A^-dependence. 

We are looking for solutions having the structure 



J AT ~ exp 



(7.36) 



where the dots stand for terms subleading in t/N. Looking at the structure of the FP 
equation, we see that all terms linear in the derivatives can be ignored at large A^ as long 
as we keep \z\ = 1. This leads us to replace the |2:| = 1 model, defined by the factors 




y, = r (7.37) 

by a new model, defined by the factors 

The Y' model preserves the U{1) symmetry of the Y model and differs from it only in the 
first-order derivative terms, while the powers of A^ appear in the same places as before. 
Therefore, the leading large- A^ behavior of the two models is the same. The advantage of 
the new model is that the Yj are restricted to an SU{1, 1) subgroup of SL{2,C), which 
forces also the product A' into SU{1, 1). This implies that the solution S'jy(A'; t) depends, 
in addition to t, on only two real variables. 

We observe now that with the right choice of variables the second-order derivatives 
acting on Sjv(A;t) only attack two of the five real arguments (on which Sjv depends, in 
addition to t) also in the Y model. Therefore, in the large- A^ limit, one can again look for a 
solution of the form (7.36). Further, this discussion indicates that significant simplifications 
will occur at large A'^ even for \z\ ^ 1, when first-order derivative terms that matter also 
in the large- A^ limit appear. 
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7.2.3 The case of \z\ = 1: The Fokker-Planck equation 

We now focus on the SU (1, 1) model and for notational convenience drop the primes, which 
now get a different use. The recursion relation determining the FP equation is 

A = A'Y (7.39) 

with 

Y = ^=== K : 1 , = £C • (7.40) 




The structure of the FP equation will be simpler in a well-chosen parametrization. The 
best choice of parameters is determined by the symmetries obeyed both by the equation 
and by our particular initial condition. 

The recursion relation has an obvious invariance A — SA, A' — SA' with S G 
SU{1, 1). This means that the FP equation for A will be invariant under multiplication 
from the left. In addition, because of the symmetry of the distribution of Y under conjuga- 
tion by the diagonal ^7(1) subgroup of SU{1, 1), there is an invariance of the equation under 
right multiplication by these group elements. Both of these invariances are broken by the 
initial condition and therefore parameters that would be associated with these invariances 
cannot be eliminated. 

However, the diagonal C/(l) subgroup of the left SU{1,1) combined with the right 
U{1) produces an U{1) x U{1) subgroup, and the initial condition is invariant under its 
diagonal subgroup. Thus, there is only one U{1) under which both the equation and the 
initial condition arc invariant. This allows us to eliminate one out of the three variables 
parametrizing SU{1, 1), making the operator O on the right-hand side of the FP equation 
a second-order partial differential operator in two real variables. However, the form of 6 is 
much more restricted, because does not know about the initial condition, and obeys many 
more symmetries. G must descend from the invariant Laplacian on the S'?7(l,l) group 
manifold after the elimination of the one extra variable. The structure of the Laplacian 
must be such that this restriction is compatible with the partial differential equation based 
on O. We do not get evolution on a coset 5J7(1, 1)/C/(1) because the variable we can 
eliminate is determined by the initial boundary condition, which singles out a specific U{1) 
subgroup that couples the right and left invariances. 

A direct derivation of the FP equation confirms the considerations above. We parametrize 
A by 

A. (;;,). |„P-|*P=1, (7.41) 

with 

a = ^e"^, b = Vu-le'''', oo>m>1, - tt < < tt . (7.42) 

In these variables the invariant measure on SU{1,1) is dudcpdip, up to a constant. The 
recursion relation for A is 



6a b 
Sb* a* 




(7.43) 
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where 5a = a — a' and 6b = b — b' . Working out the algebra, and keeping only terms up to 
second order in lo and among those only terms that could contribute to a term of the form 
|a;p, we get 

Su = -{2u - l)\ujf + ^/u{u - 1) (^e^('^-<^)cc;* + e'^^'^-'^^tc;) (7.44) 

and 

6^ = 2 [e^'^*-^^^ - e^(^-^'>u;*) . (7.45) 



The FP equation is then 



N — = -HT.n[u, (j); t) , (7.46) 



where H is a non-negative hermitian operator, 



d , ^. d u — 1 d"^ 

H = - u{u-l)- ^TTTs- 7.47 

au ou 4n am'^ 



Transforming variables to u = (1 + a;)/2, x > 1, we get 



dx^ ' dx \l{x^\) 4 

This equation is almost identical to Eq. (28) in [14]. The difference is the 1/4 term in the 
prefactor of the second derivative with respect to ^. The invariances of the equation are 
the reason for the similarity. In [14] the multiplicative random ensemble consists of real 
2x2 matrices of the form 1 + eX, where X is real and drawn from identical Gaussian 
distributions for each of its four entries. By factoring the determinant this would provide 
an evolution on the group 5L(2,]R) or, more precisely, 50(2,1). 5*^(2, M) and 5'?7(1,1) 
are the same for our continuous evolutions emanating from the identity matrix. The left 
invariances are the same in both cases since they are independent of the distribution of the 
individual factors. The difference we have obtained must therefore reflect the difference in 
the right invariances, which do depend on the distributions of the individual factors. Thus, 
one gets slightly different restrictions of the heat-kernel equation. 

The initial condition limj_^o+ Tij\i{u, (j),t) = 25{u — l)527r(<^) ((^27r('^) = 2^ Snez ^*"'^) 
is A^-in dependent, and we take the integration measure to be dudcj). Taking into account 
Eq. (7.46), we conclude that the dependence of Sjv on A'' and t is of the form 

S;v(^x,<^;t) = s(^7x,(/<;^^ . (7.49) 

Therefore, the large- limit is dominated by the short time (t = t/N) behavior of the 
probability distribution S(u, (j); t). 

In [14] the authors solve their equation by separation of variables. The (/>-dependence 
must be periodic, and is labeled by an integer m G Z. In each sector, H is replaced by 

Hr. = H = -Ux- -l)l + m^ - ]) ■ (7.50) 



dx^ 'dx \2{x + l) 4 
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The eigenfunctions and eigenvalues of Hm are known exactly. We see that our problem will 
be solved in an identical way, only the eigenvalues have to be shifted by (m/2)-^. This shift 
does not affect the matching onto the initial condition, which is the same here as in [14]. 
Therefore a formula for T,{u,(j)',t) is available, and we know that, although explicit, it is 
difficult to do much with it at the analytic level. 

7.2.4 The case of \z\ = 1: Large- AT limit from the Fokker-Planck equation 

For the modified model, Eq. (7.34) looks as follows for z = e**. 



Q{z,z*;t) 



duj #S 2^(^-cos^)^. (7.51) 



As t/N — > 0, S must become a delta-function in u and 4>. Therefore, for small t/N we 
expect S(m, 0; t/N) to drop rapidly as u increases beyond 1 and cp departs from 0. Looking 
at H, we realize that for u close to 1 the (/)-derivative term is suppressed. This leads us to 
the simple ansatz 

s(n,<A;^) ~y52.(0)e-T(«-i). (7.52) 

When this is inserted into the expression for Q the integral over (p is trivial, leaving only the 
integral over u. The latter will be dominated by a saddle point or by the endpoint u = 1. 
When the endpoint dominates we get the holomorphically factorized answer 1 1 — e** | we 
have seen before. Thus, "saddle A" corresponds to endpoint dominance. The saddle-point 
equation for u is 

- = [y/u — cos ip) ^/u , (7.53) 
and its positive solution is given by 

= i ^cos ifj + v^cos^ V' -h 2t^ . (7.54) 

We see that this saddle will first become available at a point z = e** on the unit circle only 
when 

^/u^>l, cos4'>l--, (7.55) 

in agreement with our findings earlier. Once the saddle is away from the endpoint the 
ansatz form of T, is no longer even plausible, and a more complete analysis is needed. 



7.3 The generahzed Gaussian model: Exact map to a random multiplicative 
model of 2 X 2 matrices 

When LOi ^ UJ2, that is in the generalized case, one can again reduce the problem to a 
product of random 2x2 matrices, albeit of a slightly more complicated structure than the 
one we have seen in the loi = uj2 case earlier. 
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Using similar manipulations one derives the following representation, 



/ft 



ij=i 



N 



detil-HiAj'B,] 

V j=i 



N 



(7.56) 



Here, 



and 



-UJ — UJ- 



uj-t 



(7.57) 
(7.58) 



Now, one can proceed to take the e ^ limit, deriving a FP equation for the new 2x2 
random matrix product of Aj^Bj. The structure is similar to the one in the special case 
analyzed before, and no progress can be made before the special case is fully solved. 

7.4 Large- AT universality 

The main objective of the attempt to go beyond the infinite- saddle-point approximation 
is to identify a universality class for the large- phase transition, its exponents and its 
associated relevant perturbations. For the case of (det(2; — W)) in the unitary case, this 
has been achieved in previous work. For complex matrices we need the more complicated 
object (I det(2: — M^)P), which both has a large- phase transition and a region where 
large- factorization is not useful in the sense that (| det(z — VF)P) 7^ |(det(2; — 1^))^. We 
suspect that the large- A^ universal region will have to deal with both of these issues. A 
direct approach similar to the unitary case seems difficult, but it is clear that one can make 
simplifications that do not matter at large N without loosing the universal properties. We 
have not yet learned how to do this effectively and reliably. 

A simpler case might be when ujit <^ 1. In that case we are close to the unitary case, 
with the unit circle slightly expanded into a strip of similar shape in the complex plane. 
This case should be easier to treat, in the sense of establishing large- A^ universality in an ap- 
propriately defined regime of "weak non-unit arity" . This would be an analog of the regime 
of weak non-hermiticity in non-multiplicative random complex matrix ensembles [25] . 



8. Summary 

The main objective of this paper was to generalize the universal large- A^ phase transition 
occurring in field theory models, which have important non-local observables that can be 
thought of as fluctuating unitary matrices, to the case where the non-local observables are 
more naturally chosen as complex matrices. In both cases a basic multiplicative structure 
is assumed, where the observable is a matrix given by the product of many matrices close 
to unity. The individual factors are sufficiently decorrelated to undergo a large- A^ phase 
transition in the same large- A^ universality class as simple multiplicative matrix models. 
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where the factors are completely uncorrelated. The phase transition occurs simply because 
the effective number of factors grows beyond a certain limit when a regime of strong 
coupling is entered. The effective number of factors depends both on the true number of 
factors and the departure from unity of each individual factor. 

When the coupling is small, perturbation theory indicates that parallel transport round 
a closed loop will be a matrix close to unity; when the coupling is large, the parallel 
transport is likely to depart from unity in a substantial manner. At infinite N the SL{N, C) 
random matrix associated with parallel transport round a closed finite smooth curve will 
have a spectrum which does not reach the origin of the complex plane, nor does it run away 
to infinity. Thus, we can ask if the support of the spectrum separates zero from infinity or 
not. It is possible that for small couplings it does not, but for large couplings it does. As 
the coupling is varied from zero upwards a large- phase transition would then occur. Just 
like in the unitary case, traces of finite powers of the Wilson loop matrix need not exhibit 
any nonanalyticity, but the eigenvalue distribution of the Wilson loop matrix is singular 
at the transition. 

A natural question is then what happens in a conformal gauge theory, where there 
is no difference between small and large loops and where Wilson loops corresponding to 
complex matrices naturally enter. It is quite possible that in very special cases [26] no 
large-A phase transition occurs. In such cases the behavior of the Wilson loop is captured 
by single matrix models, and there is no hint of non-commutativity playing a role. 

The conformal case can have a dimensionless coupling constant which is scale-indepen- 
dent and a free parameter of the theory. It is quite possible that less special Wilson loops 
do exhibit the above transition. 

Even in the unitary cases it may be advantageous to use a regularization that introduces 
some weak non-unit arity, in which case a study of the universality of the appropriate matrix 
model would be useful. 

While our work shows beyond reasonable doubt that the Durhuus-Olesen transition 
of the unitary matrix case generalizes to complex matrices, more work is needed to bring 
the understanding of large-A universality in the complex case to the same level as that of 
the unitary case. While the latter corresponds to a regularization [9] of the singularity- 
generating mechanism of the inviscid Burgers equation for wave propagation restricted to 
some line, we now need something similar for the inviscid Burgers equation for complex 
waves propagating in the complex plane. Much is known about this, and it is hard to believe 
that the particular regularization employed by going to the large- A^ "double-scaling" limit 
will generate something totally new. To know what exactly happens more work is needed, 
but the vast existing applied mathematics literature indicates that the problem is soluble. 
There exists an example in the context of single matrix models where a connection to the 
deformation theory of some specific partial differential equations has been established [27]. 
We plan to return to this topic in the future. 
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